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Abstract 

The kinetic theory is formulated with respect to anholonomic frames of 
reference on curved spacetimes. By using the concept of nonlinear connection 
we develop an approach to modelling locally anisotropic kinetic processes and, 
in corresponding limits, the relativistic non-equilibrium thermodynamics with 
local anisotropy. This lead to a unified formulation of the kinetic equations 
on (pseudo) Riemannian spaces and in various higher dimensional models of 
Kaluza-Klein type and/or generalized Lagrange and Finsler spaces. The tran- 
sition rate considered for the locally anisotropic transport equations is related 
to the differential cross section and spacetime parameters of anisotropy. The 
equations of states for pressure and energy in locally anisotropic thermody- 
namics are derived. The obtained general expressions for heat conductivity, 
shear and volume viscosity coefficients are applied to determine the transport 
coefficients of cosmic fluids in spacetimes with generic local anisotropy. We 
emphasize that such locally anisotropic structures are induced also in gen- 
eral relativity if we are modelling physical processes with respect to frames 
with mixed sets of holonomic and anholonomic basis vectors which naturally 
admits an associated nonlinear connection structure. 
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1 Introduction 



The experimental data on anisotropies in the microwave background radiation (see, 
for instance, Ref. and [0) and modern physical theories support the idea 



that in the very beginning the Universe has to be described as an anisotropic and 
higher dimension spacetime. It is of interest the investigation of higher dimension 
generalized Kaluza-Klein and string anisotropic cosmologies. Therefore, in order to 
understand the initial dynamical behavior of an anisotropic universe, in particular, to 
study possible mechanisms of anisotropic inflation connected with higher dimensions 
[0 we have to know how we can compute the parameters (transport coefficients, 
damping terms and viscosity coefficients) that characterize the cosmological fluids 
in spacetimes with generic anisotropy. 

The first relativistic macroscopic thermodynamic theories have been proposed 
in Refs. H and p2|. The further developments and applications in gravitational 



physics, astrophysics and cosmology are connected with papers |2^, ^ |24 



The Israel's approach to a microscopic Boltzmann like kinetic theory for relativistic 



gases |]T8| makes possible to express the transport coefficients via the differential 
cross sections of the fluid's particles. Here we also note the Chernikov's results on 
Boltzmann equations with collision integral on (pseudo) Riemannian spaces 0. A 
complete formulation of relativistic kinetics is contained in the monograph |]l5l by 
de Groot, van Leeuwen and van Weert. We also emphasize the Vlasov's monograph 
3^ were an attempt to statistical motivation of kinetic and thermodynamic theory 
on phase spaces enabled with Finsler like metrics and connection structures was 
proposed. Recent developments and applications of kinetics and nonequilibrium 
thermodynamics could be found in PS[ . 

The extension of the four dimensional considerations to higher dimensions is due 
to 1 29 1 and [Q. The generalization of kinetic and thermodynamic equations and 
formulas to curved spaces is not a trivial task. In order to consider flows of par- 
ticles with noninteger spins, interactions with gauge flelds, and various anisotropic 
processes of microscopic or macroscopic nature it is necessary a reformulation of the 
kinetic theory in curved spacetimes by using the Cartan's moving frame method 
[^]. A general approach with respect to anholonomic frames contains the possi- 
bility to take into account generic spacetime anisotropies which play an important 
role in the vicinity of cosmological or astrophysical singularities and for non-trivial 
reductions of higher dimension theories to lower dimensional ones. In our works 
we proposed to describe such locally anisotropic spacetimes and inter- 
actions by applying the concept of nonlinear connection (in brief, N-connection) 
fleld which models the local splitting of spacetime into horizontal (isotropic) and 
vertical (anisotropic) subspaces. For some values of components the N-connection 
could parametrize, for instance, toroidal compactiflcations of higher dimensions but, 
in general, its dynamics is to be determined, in different approaches, by fleld equa- 
tions or constraints in some generalized Kaluza-Klein, gauge gravity or (super)string 
theories 130,131102153. 



The concept of N-connection was flrstly applied in the framework of Finsler and 
Lagrange geometry and gravity and their higher order extensions (see Refs. p3| , pT] , 
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|5D| , on details). Here it is to be emphasized that every generahzed Finsler hke 
geometry could be modelled equivalently on higher dimensions (pseudo)Riemannian 
spacetimes (for some models by introducing additional torsion and/or nonmetricity 
structures) with correspondingly adapted anholonomic frame structures. If we re- 
strict our considerations only to the higher dimensional Einstein gravity, the induced 
N-connection structure becomes a "pure" anholonomic frame effect which points to 
the fact that the set of dynamical gravitational field variables given by metric's com- 
ponents was redefined by introducing local frame variables. The components of a 
fixed basis define a system (equivalently, frame) of reference (in four dimensions one 
uses, for instance, the terms of vierbien, or tetradic field) with respect to which, in 
its turn, one states the coefficients of curved spacetime's metric and of fundamental 
physical values and field equations. It should be noted that the procedure of choos- 
ing (establishing) of a system of reference must be also physically motivated and 
that in general relativity this task is not considered as a dynamical one following 
from the field equations. For trivial models of physical interactions on curved space- 
times one can restrict our considerations only with holonomic frames which locally 
are linearly equivalent to some coordinate basis. Extending the class of physical 
fields and interactions (even in the framework of Einstein's gravity), for example, 
by introducing spinor fields, statistical and fluid models with spinning particles we 
have to apply frame bundle methods and deal with general anholonomic frames and 
phase spaces provided with Cartan's Finsler like connections (induced by the Levi 
Civita connection on (pseudo)Riemannian spacetime) 0. 

The modelling of kinetic processes with respect to anholonomic frames (which 
induces corresponding N-connection structures) is very useful with the aim to elu- 
cidate flows of fluids of particles not being in local equilibrium. More exactly, such 
generic anisotropic fluids are considered to be like in a thermodynamic equilibrium 
with respect to some frames which are locally adapted both to the spacetime met- 
ric and N-connection structures but in general (for instance, with respect to local 
coordinate frames) the conditions of local equilibrium are not satisfied. 

We developed a proper concept of locally anisotropic spacetime (in brief, la- 
spacetime) [^, ^ |3^ in order to provide a unified (super) geometric background 
for generalized Kaluza-Klein (super) gravities and low energy (super) string models 
when the higher dimension spacetime is characterized by generic local anisotropics 
and compatible metric, N-connection and correspondingly adapted linear connec- 
tion structures. It should be noted that the term la-spacetime can be used even 
for Einstein spaces if they are provided with anholonomic frame structures. Our 
treatment of local anisotropy is more general than that from Refs. which is used 
for a subclass of Finsler like metrics being conformally equivalent (with conformal 
factors depending both on spacetime coordinates and tangent vectors) to the flat 
(pseudo) Euclidean, equivalently, Minkovschi metric. If in the Bogoslovsky, Goenner 
and Asanov works |lT| on Finsler like spacetimes there are investigated possible 
effects of violation of local Lorentz and/or Poincare symmetries, our approach to 
locally anisotropic strings, field interactions and stochastics ^ |3^ is back- 
grounded on the fact that such models could be constructed as to be locally Lorentz 
invariant with respect to frames locally adapted to N-connection structures. 
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The purpose of the present article is to formulate a Boltzmann type kinetic 
theory and non-equilibrium thermodynamics on spacetimes of higher dimension 
with respect to anholonomic frames which models local anisotropics in both type 
Einstein or generalized Finsler-Kaluza-Klein theories. We shall also discuss possible 
applications in modern cosmology and astrophysics. 

On kinetic and thermodynamic part our notations and derivations are often in- 
spired by the Refs. 0] and |jl5| . A number of formulas will be similar for both locally 
isotropic and anistropic spacetimes if in the last case we shall consider the equa- 
tions and vector, tensor, spinor and connection objects with respect to anholonomic 
rest frames locally adapted to the N-connection structures. In consequence, some 
tedious proofs and intermediary formulas will be omitted by referring the reader to 
the corresponding works. 

To keep the article self-consistent we start up in Section 2 with an overview 
of the so-called locally anisotropic spacetime geometry and gravity. In Sections 3 
and 4 we present the basic definitions for locally anisotropic distribution functions, 
particle flow and energy momentum tensors and generalize the kinetic equations for 
la~spacetimes. The equilibrium state and derivation of expressions for the particle 
density and the entropy density for systems that are closed to a locally anisotropic 
state of chemical equilibrium are considered in Section 5. Section 6 is devoted to the 
linearized locally anisotropic transport theory. There we discuss the problem of solu- 
tion of kinetic equations, prove the linear lows for locally anisotropic non-equilibrium 
thermodynamics and obtain the explicit formulas for transport coefficients. As an 
example, in Section 7, we examine the transport theory in curved spaces with rota- 
tion ellipsoidal horizons. Concluding remarks are contained in Section 8. 



2 Spacetimes with Local Anisotropy 

We outline the necessary background on anholonomic frames and nonlinear con- 
nections (N-connection) modelling local anisotropics (la) in curved spaces and on 
locally anisotropic gravity |Q |31[] (see Refs. and for details on space- 

time differential geometry and N-connections structures). We shall prove that the 
Cartan's moving frame method allows a geometric treatment of both type of 
locally isotropic (for simplicity, we shall consider (pseudo) Riemannian spaces) and 
anisotropic (the so-called generalized Finsler-Kaluza-Klein spaces). 



2.1 Anholonomic frames and Einstein equations 

In this paper spacetimes are modelled as smooth (i.e class C°°) manifolds V(^d) of finite 
integer dimension d > 3,4,..., being Hausdorff, paracompact and connected. We 
denote the local coordinates on V(^d) by variables u", where Greek indices ... = 
3,4, ... could be both type coordinate or abstract (Penrose's) ones. A spacetime is 
provided with corresponding geometric structures of symmetric metric gap and of 
linear, in general nonsymmetric, connection F"^^ defining the covariant derivation 
Va satisfying the metricity conditions V aQp-y = 0. We shall underline indices, a, /3, 
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if one would like to emphasize them as abstract ones. 

Let a set of basis (frame) vectors e- = {e§ = S'o/je-'^} on V(rf) be numbered by 
an underlined index. We shall consider only frames associated to symmetric metric 
structures via relations of type 

e^e^a = rja,i3_ and e^e^^£ = 9ai3, 

where the Einstein summation rule is accepted and rja^p is a given constant symmetric 
matrix, for simplicity a pseudo-Euclidean metric of signature (— , +....+) (the sign 
minus is used in this work for the time like coordinate of spacetime). Operations 
with underlined and non-underlined indices are correspondingly performed by using 
the matrix r]a,/3, its inverse r]^, and the metric ga/s and its inverse g""^. A frame 
(local basis) structure Ca on V(^d) is characterized by its anholonomy coefficients w"^^ 
defined from relations 

eaC/j - C/jCa = w\pe^. (1) 

With respect to a fixed basis and its dual we can decompose tensors and 
write down their components, for instance, 

where by ® it is denoted the tensor product. 

A spacetime V(^d) is holonomic (locally integrable) if it admits a frame 
structure for which the anholonomy coefficients from (|l]) vanishes, i.e. w'^^^ = 0. In 
this case we can introduce local coordinate bases, 

= d/du'' (2) 

and their duals 

= du"" (3) 

and consider components of geometrical objects with respect to such frames. 

We note that the general relativity theory was formally defined on holonomic 
pseudo-Riemannian manifolds. Even on holonomic spacetimes, for various (geo- 
metrical, computational and physically motivated) purposes, it is convenient to use 
anholonomic frames e^, but we emphasize that for such spacetimes one can always 
define some linear transforms of frames to a coordinate basis, = CL^da'- By 
applying both holonomic and anholonomic frames and theirs mutual transforms 
on holonomic pseudo-Riemannian spaces there were developed different variants of 
tetradic and spinor gravity and extensions to linear, affine and de Sitter gauge group 
gravity models [|1^, [1^. 

A spacetime is generically anholonomic (locally non-integrable) if it does 
not admit a frame structure for which the anholonomy coefficients from (|I]) vanishes, 
i.e. w^Q,^ 7^ 0. In this case the anholonomy becomes a proper spacetime characteris- 
tics. For instance, a generic anholonomy could be obtained if we consider nontrivial 
reductions from higher dimension spaces to lower dimension ones. It induces non- 
vanishing additional terms into the torsion, 

T {6^,6/3) = T"'^ Ja, 
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and curvature, 

tensors of a linear connection T"^^, with coefficients defined respectively as 

T"^^ = r V - r V + (4) 

and 

-^/s^T = ^-r^^p-y ~ ^i^'^i35 + r'^^^r"^^^ — T'^p^T"'^^ + r°g^w'^^^. (5) 

The Ricci tensor is defined 

-R/37 = Rft-ya (6) 



and the scalar curvature is 



R = g^^Rp,. (7) 



The Einstein equations on an anholonomic spacetime are introduced in a stan- 
dard manner, 

Rp-r ~ 1^9 mR = kTp^, (8) 

where the energy-momentum d-tensor T^^ includes the cosmological constant terms 
and possible contributions of torsion (H) and matter and k is the coupling constant. 
For a symmetric linear connection the torsion field can be considered as induced by 
some anholonomy (or equivalently, by some imposed constraints) conditions. For 
dynamical torsions there are necessary additional field equations, see, for instance. 



the case of locally anisotropic gauge like theories |^ . 

The usual locally isotropic Einstein gravity is obtained on the supposition that 
for every anholonomic frame could be defined corresponding linear transforms to a 
coordinate frame basis. 

It is a topic of further theoretical and experimental investigations to establish if 
the present day experimental data on anisotropic structure of Universe is a conse- 
quence of matter and quantum fluctuation induced anisotropies and for some scales 
the anisotropy is a consequence of anholonomy of observer's frame. The spacetime 
anisotropy could be also a generic property following, for instance, from string the- 
ory, and from a more general self-consistent gravitational theory when both the left 
(geometric) and right (matter energy-momentum tensor) parts of Einstein equations 
depend on anisotropic parameters. 



2.2 The local anisotropy and nonlinear connection 

A subclass of anholonomic spacetimes consists from those with local anisotropy 
modelled by a nonlinear connection structure. In this subsection we briefly outline 
the geometry of anholonomic frames with induced nonlinear connection structure. 

The la-spacetime dimension is split locally into two components, n for isotropic 
coordinates and m for anisotropic coordinates, when n(^a) = n + m with n > 2 and 
m > 1. We shall use local coordinates = (x*,?/"), where Greek indices ... 
take values 1,2, ...,n + m and Latin indices i and a are correspondingly n and m 
dimensional, i.e. k... = 1, 2, n and a, b, c, ... = 1, 2, m. 
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Now, we consider an invariant geometric definition of spacetime's splitting into 
isotropic and anisotropic components. For modelling la-spacetimes one uses a vec- 
tor bundle E = {En+m,P, M(^n), F(^m),Gr) provided with nonlinear connection (in 
brief, N— connection) structure = {Nj where Nj (m") are its coefficients 



2^ . We use denotations: Em+n is a (n + m)-dimensional total space of a vector 
bundle, M(„) is the ra-dimensional base manifold, is the typical fiber being a 
m-dimensional real vector space, Gr is the group of automorphisms of and p is 
a surjective map. For simplicity, we shall consider only local constructions on vector 
bundles. 

The N-connection is a new geometric object which generalizes that of linear 
connection. This concept came from Finsler geometry (see the Cartan's monograph 
[§]), the global formulation of it is due to W. Barthel 0, and it is studied in details 



in Miron and Anastasiei works p3[. We have extended the geometric constructions 



for spinor bundles and superbundles with further applications in locally anisotropic 
field theory and strings and modern cosmology and astrophysics pOl, pll, [33]. We 



have also illustrated that the N-connection could be introduced on (pseudo) 
Riemnannian spacetimes and in Einstein gravity if we consider anholonomic frames 
consisting from sets of basic vectors some of them being holonomic and the rest 
anholonomic. In this case the N-connection coefficients are associated to the frame 
structure and transforms into some metric components if the considerations are 
transferred with respect to a coordinate basis. 

The rigorous mathematical definition of N-connection is based on the formalism 
of horizontal and vertical subbundles and on exact sequences of vector bundles. 
Here, for simplicity, we define a N-connection as a distribution which for every 
point u = {x,y) G £ defines a local decomposition of the tangent space of our vector 
bundle, TuE, into horizontal, HuE, and vertical (anisotropy), VuE, subspaces, i.e. 

T^E = HuE © VuE. 

If a N-connection with coefficients Nj {u'^) is introduced on the vector bundle 
S the modelled spacetime posses a generic local anisotropy and in this case we can 
not apply in a usual manner the operators of partial derivatives and their duals, 
differentials. Instead of coordinate bases (0) and (H) we must consider some bases 
adapted to the N-connection structure: 

= (5„ da) = ^= {6^ = ^ = 4- -Nl: (x^, y") St, da = i^] (9) 

and 

= {d\ 5") = 6u^ = (ct = dx\ 5" = 6y^ = dy" + (x\y^) dx^) . (10) 
A nonlinear connection (N-connection) is characterized by its curvature 

% = 4^--^ + ^f^T - ^f^T- (11) 
dx^ dx' ' dy^ ^ dy^ ^ ^ 
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Here we note that the class of usual linear connections can be considered as a 
particular case when 

The elongation (by N-connection) of partial derivatives in the adapted to the N- 
connection partial derivatives (^), or the locally adapted basis (la-basis) reflects 
the fact that the spacetime £ is locally anisotropic and generically anholonomic 
because there are satisfied anholonomy relations (|I]), 

where anolonomy coefficients are as follows 

= 0,w\j = 0,w\, = 0,w'',, = 0,w',, = 0, 

On a la-spacetime the geometrical objects have a distinguished (by N-connec- 
tion), into horizontal and vertical components, character. They are briefly called 
d-tensors, d-metrics and/or d-connections. Their components are defined with 
respect to a la-basis of type (H), its dual (0), or their tensor products (d-linear 
or d-affine transforms of such frames could also be considered). For instance a 
covariant and contravariant d-tensor Z, is expressed as 

Z = Z%5^ = Z'fii ® & + Z;\bi ® 5'^ + Z\f)y, ® + Z\db ® 5\ (12) 

A symmetric d-metric on la-space £ is written as 

^s^ = gap {u) r ® 5^ = gij{x, y)dx'dx^ + hab{x, y)dy''Sy\ (13) 

A linear d-connection D on la-space £, 

D5,5/3 = 0^ {x',y)s^, 

is parametrized by non-trivial h-v-components, 

^ ^*ic5 ^""bc) ■ (14) 

Some d-connection and d-metric structures are compatible if there are satisfied the 
conditions 

Dag/Sy = 0. 

For instance, a canonical compatible d-connection 

c-pa I'^T^ era c /^i c r~ia \ 

^ /37 ~ jfcJ fefc' ^ jc^ ^ be) 

is defined by the coefficients of d-metric (p!3D, gij {x,y) and hab {x,y) , and by the 
coefficients of N-connection, 



e T i 
^ jk - 


^g'"" i^kgnj + Sjgnk ■ 






era 
^ bk — 




- hdcdbNf - 


- hdbdcNf) , 


c/~ii 

jc - 


\g'^dcgjk, 






e/^a 

be — 


hi""" {d.hdb + dbhdc 


- ddhbc) 
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J-.jk 


rpl 

- J-jk - 


rjii 
•ja 


= 0, T 


rpa 

•ij 


— -f]". 



The coefficients of the canonical d-connection generahze for la-spacetimes the well 
known Cristoffel symbols. 

For a d-connection (|T^) we can compute the components of, in our case d-torsion, 



jk kji -'-ja ^.jai-'-aj ja' 

be — '-'.be — '-'be '-"cfe; 

rpa o Ara r a rpa rpa 

-I- M — '^bJ-^ i ~^.bj-> Ab — M- 

In a similar manner, putting non-vanishing coefficients into the formula for 
curvature (|^), we can compute the non-trivial components of a d-curvature 

r}.i )\ T ^ Ar^ J- T ^ T ^ jm ji /^i c\a 

^h.jk — '^kJ^.hj 'Jj-'^.hk ~^ -'^.hj-'^mk ^.hk^mj ^.ha^'-.jk^ 

^b.jk = ^kL%j — + L^i^jL'^ck ~ ^bk^^ ~ C°bc^''.jki 

Pf.ka = ^kL^.jk + Cji,T^i^^ — (9fcCj„ + Uj^C^jg^ — L\^f^C\^ — L'^,ak^^.jc}-i 

jD-e Pi T c I r^e rpd r-ic \ TC ryd jd r^c jd r-^e \ 

^b.ka — ^ai-^.bk + .bd^ .ka ~ \^k^ .ba + ^ .dk^ .ba ^ ^ .bk^ .da ~ ^ .ak^ .bd) 

Q.i f\ /^i a /^i I /^h /^i /^h /^i 

•^j.bc ~ ^c^.jb .je^ ^ .jb^ .he ^.jc^hbi 

Q.a o r-ia o /~ia < /~ie /^a /~ie r-ia 

^b.ed — ^d^.bc ~ (^e^M + ^ .bc^ .ed ~ ^ .bd^ .ec' 



The components of the Ricci tensor (|T6D with respect to locally adapted frames 



(H) and (^) (in this case, d-tensor) are as follows: 

Rij = Ri,jki Ria = ~ Pia = ^Pi.kai (IS) 

p 1 p p.b p q.c 

J^ai — J^ai — a.ib^ ~ '-'a.bc- 

We point out that because, in general, ^Pai 7^ the Ricci d-tensor is non sym- 
metric. 

Having defined a d-metric of type ([T3|) in S we can compute the scalar curvature 
(|^ of a d-connection D, 

= C'^R^p = R + S, (17) 
where R = g'^ij and S = h^^Sab- 



Now, by introducing the values (|TB]) and (0) into anholonomic gravity field 
equations (§) we can write down the system of Einstein equations for la-gravity 
with prescribed N-connection structure 



Ri,-]^(R + S)g,, = kT,^, (18) 

Sab ~ 2 ("^ ~'~ ^"'^ ^ ^^ab) 

1 p . _ . 



2 



Pia i 



where Tj^, Ta^, T^j and Tj^ are the components of the energy-momentum d-tensor 
field. We note that such decompositions into h- and v-components of gravitational 
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field equations have to be considered even in general relativity if physical interactions 
are examined with respect to an anholonomic frame of reference with associated N- 
connection structure. 

There are variants of la-gravitational field equations derived in the low-energy 
limits of the theory of locally anisotropic (super) strings |3ll] or in the framework 
of gauge like la-gravity |^5|, ^ when the N-connection and torsions are dynamical 
fields and satisfy some additional field equations. 



2.3 Modelling of generalized Finsler geometries in 
(pseudo) Riemannian spaces 

The present day trend is to consider the Finsler like geometries and their general- 
izations as to be quite sophisticate for straightforward applications in quantum and 
classical field theory. The aim of this subsection is to proof that, a matter of prin- 
ciple, such geometries could be equivalently modelled on corresponding (pseudo) 
Riemannian manifolds (tangent or vector bundles) by using the Cartan's moving 
frame method and from this viewpoint a wide class of Finsler like metrics could be 
treated as some solutions of usual Einstein field equations. 



2.3.1 Almost Hermitian Models of Lagrange and Finsler Spaces 

This topic was originally investigated by Miron and Anastasiei ; here we outline 
some basic results. 

Let us model a la-spacetime not on a vector bundle £ but on a manifold TM = 
TM\{0} associated to a tangent bundle TM of a n-dimensional base space M (when 
the dimensions of the typical fiber and base are equal, n = m and \{0} means that 
there is eliminated the null cross-section of the bundle projection r : TM M) 
and consider d-metrics of type 

Ss"^ = ga(3 (u) 5"0 5^ = gij{x, y)dx'dx^ + gij{x, y)Sy'5y\ (19) 
On TM we can define a natural almost complex structure C(a) as follows 

C(a) {S^) = -d/dy' and q,) [d/dy') = 5, 

where the la-derivative 5i = d/dx^ — Nfd/dy^ (|^) and la-differential 5* = dy^ + 
N^dx^ (jiy) act on TM being adapted to a nontrivial N-connection structure N = 
{Nj {x,y)} in TM. It is obvious that C'^^-^ = —I. The pair (^5s^,C(a)) defines an 
almost Hermitian structure on TM with an associate 2-form 

e = (x, y) 5'Adx^ 

and the triad K = (tm, 

6s'^,Cta)) is an almost Kahlerian space. By straightfor- 
ward calculations we can verify that the canonical d-connection (|T^) satisfies the 
conditions 



10 



for any d-vector X on TM and has zero hhh- and f f f-torsions. 



The notion of Lagrange space [EH, |23] was introduced as a generahzation of 



Finsler geometry in order to geometrize the fundamental concepts in mechanics. A 
regular Lagrangian L (x*, y*) on TM is introduced as a continuity class C°° function 
L : TM IR for which the matrix 

g (x,y) = -—— (20) 

has rank n and is of constant signature on TM. A d-metric ([l9|) with coefficients of 
form (pOD, a corresponding canonical d-connection ( p!^ and almost complex struc- 
ture C(a) defines an almost Hermitian model of Lagrange geometry. 

For arbitrary metrics g-ij (x, y) of rank n and constant signature on TM, which 
can not be determined as a second derivative of a Lagrangian, one defines the so- 



called generalized Lagrange geometry on TM (see details in [^). 

A particular subclass of metrics of type (^) consists from those where instead 
of a regular Lagrangian one considers a Finsler metric function F on M defined 
as F : TM IR having the properties that it is of class C°° on TM and only 
cjDntinuous on the image of the null cross-section in TM, the restriction of F on 
TM is a positive function homogeneous of degree 1 with respect to the variables 
i. e. 

T(x,Ay) = AT (x,y),AGlR^ 
and the quadratic form on H" with coefficients 

1 a^T^ 

defined on TM, is positive definite. Different approaches to Finsler geometry, its 
generalizations and applications are examined in a number of monographs fll], |^, ^ 



and as a rule they are based on the assertion that in this type of geometries the 
usual (pseudo)Riemannian metric interval 



ds = \J Qij (x) dx^dx^ 

on a manifold M is changed into a nonlinear one defined by the Finsler metric 
F (fundamental function) on TM (we note an ambiguity in terminology used in 
monographs on Finsler geometry and on gravity theories with respect to such terms 
as Minkowschi space, metric function and so on) 

ds = F (^x\ dx^) . (22) 

Geometric spaces with a 'combersome' variational calculus and a number of curva- 
tures, torsions and invariants connected with nonlinear metric intervals of type (pO|) 
are considered as less suitable for purposes of modern field and particle physics. 
In our investigations of generalized Finsler geometries in (super) string, gravity 



and gauge theories |31| we advocated the idea that instead of usual geometric 
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constructions based on straightforward applications of derivatives of (^) following 
from a nonlinear interval (|22| ) one should consider d-metrics (0) and/or (|1^) with 
the coefficients of necessity determined via an almost Hermitian model of a Lagrange 
(pOD, Finsler geometry ( pT] ) and/or their extended variants. This way, by a synthesis 



of the moving frame method with the geometry of N-connection, we can investigate 
a various class of higher and lower dimension gravitational models with generic or 
induced anisotropies in a unified manner on some anholonomic and / or Kaluza-Klein 
spacetimes. 

As a matter of principle, having a physical model with a d-metric and geometrical 
objects associated to la-frames, we can redefine the physical values with respect to 
a local coordinate base on a (pseudo) Riemannian space. The coefficients Qij {x, y) 
and hab {x,y) of d-metric written for a la-basis 6u°' = {dx^,6y'^) transforms into a 
usual (pseudo) Riemannian metric if we rearrange the components with respect to 
a local coordinate basis du'^ = (^dx*, dy'^^, 



^ij J j ab i I 1^23) 



be 



where 'hats' on indices emphasize that coefficients of the metric are given with re- 
spect to a coordinate (holonomic) basis on a spacetime V^"+'". Parametrizations 
(ansatzs) of metrics of type (^) are largely applied in Kaluza-Klein gravity and its 
generalizations p5| . In our works |3^, Q , following the geometric construc- 



tions from |2^, we proved that the physical model of interactions is substantially 
simplified, as well we can correctly elucidate anisotropic effects, if we work with 
diagonal blocks of d-metrics (|13D with respect to anholonomic frames determined 
by a N-connection structure. 

2.3.2 Finsler like metrics in Einstein's gravity 



There are obtained (see Appendix) some classes of locally anisotropic cos- 

mological and black hole like solutions (in three, four and higher dimensions) which 
can be treated as generalized Finsler metrics being of nonspheric symmetry (with 
rotation ellipsoid, torus and cylindrical event horizons, or with elliptical oscillations 
of horizons). Under corresponding conditions such metrics could be solutions of field 
equations in general relativity or its lower or higher dimension variants. Here we 
shall formulate the general criteria when a Finsler like metric could be a solution of 
gravitational field equations in Einstein gravity. 

Let consider on TM an ansatz of type (|2^) when gij = hij = ^d"^ F"^ / dy^ dy^ (for 
simplicity, we omit 'hats' on indices) i.e. 

n — _ dy^dyi i ^^jdyi'dy' ^^jdy'^dy^ \ (nA\ 

\ i dy^dy^ dy^dy^ 



A metric (p^, induced by a Finsler quadratic form ( pi]) could be treated in a frame- 
work of a Kaluza-Klein model if for some values of Finsler metric F {x, y) and N- 
connection coefficients (x, y) this metric is a solution of the Einstein equations 
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written with respect to a holonomic frame. For the dimension = 2, when the 
values F and are chosen to induce locally a (pseudo) Riemannian metric g-^-j^ of 
signature (— , +, +, +) and with coefficients satisfying the four dimensional Einstein 
equations we define a subclass of Finsler metrics in the framework of general rela- 
tivity. Here we note that, in general, a N-connection, on a Finsler space, subjected 
to the condition that the induced (pseudo) Riemannian metric is a solution of usual 
Einstein equations does not coincide with the well known Cartan's N-connection 



11]. We have to examine possible compatible deformations of N-connection 



structures PB| . 

Instead of Finsler like quadratic forms we can consider ansatzs of type (p3|) 
with Qij and hij induced by a Lagrange quadratic form (|20|) . A general approach 
to the geometry of spacetimes with generic local anisotropy can be developed on 
embeddings into corresponding Kaluza-Klein theories and adequate modelling of 
la-interactions with respect to anholonomic or holonomic frames and associated 
N-connection structures. 



3 Collisionless relativistic kinetic equation 

As argued in Section 2, the spacetimes could be of generic local anisotropy after non- 
trivial reductions from some higher dimension theories or posses a local anisotropy 
induced by anholonomic frame structures even we restrict our considerations to the 
general relativity theory. In this line of particular interest is the formulation of rel- 
ativistic kinetic theory with respect to general anholonomic frames and elucidation 
of locally anisotropic kinetic and thermodynamic processes. 



3.1 The distribution function and its moments 

We use the relativistic approach to kinetic theory (we refer readers to monographs 



1^, |3^ for hystory and complete treatment). Let us consider a simple system 
consisting from w point particles of mass m in a la-spacetime with a d-metric 
Qafi- Every particle is characterized by its coordinates m";) = (^x^;-,, , {v} = = 
ct is considered the time like coordinate, for simplicity we put hereafter the light 
velocity c = 1); x^^^, x^^^, x";-, and ?/(\), y^;-,, respectively space like and 
anisotropy coordinates, where the index (/) enumerates the particles in the system. 
The particles momenta are denoted by Pa{i) = gaisP^i), ■ We shall use the distribution 
function <l> {u^'.pp) , given on the space of supporting elements (m",^/?) , as a general 
characteristic of particle system. 

The system (of particles) is to be defined by using the random function 

</)(n",p^)=5:/535("(«))(n"-n^,)(s)) 5^^^^^^ [p^ - p^^is)) , (25) 
1=1 •' 

where the sum is taken on all system's particles, 

5s = \J gap5u'^5ui^ 
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is the interval element along the particle trajectory parametrized by a natural 

parameter s and (J^^") (-u") is the n(a)-dimensional delta function. The functions 
and describing the propagation of the /-particle are found from the 

motion equations on la-spacetime 

m ^ = vl) (26) 



ds 

ds~ ~ ~S ^ ap[U^i^[s)jPr{l)P(^i)- t^ii), 



where '^T"^^ is the canonical d-connection with coefficients (|I5|) and by Fa(i) is 
denoted an exterior force (electromagnetic or another type) acting on the /th particle. 

The distribution function $ (m,p) is defined by averaging on paths of random 
function 

$ {U,p) =« {U,p) » 

where brackets << ... >> denote path averaging. 

Let consider a space like hypersurface J-'{u") = const with elements = UaST., 

where = and dE = ^fjl'" , \Dr\ = ^g^PD^DpT and is the 

invariant space-time volume. A local system of reference in a point u^q-,, with the 

metric gf^^ = rjij = diag {—1,1,1, , is obtained if JF(M"g^) = xj^-^ . In this 
case SEa = ^^f")"^!* = dx'^dx^...dx"'Sy^Sy'^...6y'^ and 6^ is the Kronecker 

symbol. The value 

<l>(M,p)i;"dS, ^ 



9\ 

where = p°'/m and 6P = S"''-'^'>p/ ^\g\ = Spi6p2---Sprn^^/ \f\g\ is the invariant 
volume in the momentum space and defines the quantity of particles intersecting 
the hypersurface element (5Eq with momenta Pa included into the element 6P in 
the vicinity of the point u". The first < > and second moments < p'^p^ > 
of distribution function $ [x, p) give respectively the fiux of particles n" and the 
energy-momentum T"^, 

< p° >= y $ (u, p) p''6P = mn" (27) 

and 

< p"/ >= / $ {u, p) pyf^dP = mT"^. (28) 

We emphasize that the motion equations (p6| ) have the first integral, g°'^p^^pf' = 
m? = const for every /-particle, so the functions (pSf ) and, in consequence, $ {u,p) 
are non-zero only on the mass hypersurface 

= m', (29) 

which in distinguished by a N-connection form (see the dual to d-metric (pTSD) is 
written 

g'^PiPj +g''^PaPh = m^- 
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For computations it is convenient to use a new distribution function / {u,p) on 
a ?T,(a) — 1 dimensional la-hyperspace (|29| ) 

$ {u^,pp) = f {u-,p^) 6 (v^rtw^ - ^) (pt) 

where 



1, Pt>0 
0, Pr<0 



and the index (3 runs values in n(a) — 1 dimensional la-space. The flux of particles 
( p7D and of energy-momentum (^) are computed by using f{u,p) as 



with respective base and fiber n"' components, 

= 6q p'^f {vC^p^ and 72"= j 6q p" f {v? ,p^ 



when 

T"/5 = j p''p^f{u\p^) 

where p^ is expressed via ...jp^^") by using the equation (^U^ and it is used 

the abbreviation S<; = S"'^'^^~^p/{^J\g\p~). The energy-momentum can be also split 
into base-fiber (horizontal-vertical, in brief, h- and v-, or hv-components) by a 
corresponding distinguishing of momenta, T"^ = {T*-', T'^-', T*^, T'^''}. 

The one-particle distribution function / {u,p) characterizes the number of parti- 
cles (with mass m, or massless if m = 0) at a point of a la-spacetime of dimension 
n(a) having distinguished by N-connection momentum vector (in brief, momentum 

d-vector) p° = (p*,p") = {p^,J^) ■ One states that 



/ (m, p) p^ 6ap 6q = f (m, p) p^ 6ai 6q + f {u, p) p" 6aa Sq, 

gives the number of world-lines of particles with momentum d-vector p" in an 
interval Sp"' around p, crossing a space like hypersurface 5(j/3 at a point u. When 5(j/3 
is taken to be time like one has 6af3 = (^"('''m, 0, 0) . 

The velocity d-field of a fluid If^ in la-spacetime can be defined in some way 
like in relativistic kinetic theory [l^. To obtain a particularly simple form for 
the energy-momentum d-tensor one should follow the Landau-Lifshitz approach 



2^ when the locally anisotropic fluid velocity (u) = {U'^ (u) (u)) and the 
energy density e (u) are defined respectively as the eigenvector and eigenvalue of the 
eigenvalue equation 

T"^f/^ = eU". (30) 

A unique value of f/" is to be found from the conditions to be a time like and 
normalized to the unity d-vector, If^Ua = 1- 
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Having so fixed the fluid d-velocity we can define correspondingly the particle 
density 

n{u)=n'^Ua, (31) 

the energy density 

eiu) = U^T^^Up (32) 
and the average energy per particle 

e = e/n. (33) 
We can introduce also the pressure d-tensor 

{u) = A^^TruA"^ (34) 

by applying the projector 

^af} ^ gap _ ^a^p ^35^) 

with properties 

A°^A^^ = A°, A^I^Up = and A{: = n(„) - 1. 

With respect to a la-frame (|^) we can introduce a particular Lorentz system, 
called the locally anisotropic rest frame of the fluid when = (1, 0, 0, 0) and 
A° = diag (0, 1, 1, 1) . So the pressure d-tensor was defined as to coincide with the 
la-space part of the energy-momentum d-tensor with respect to a locally anisotropic 
rest frame. 



3.2 Collisionless kinetic equation 

We start our proof by applying the identity 

- «f,)(s)) (/ - pf,)(.)) = 0. 

Differentiating under this integral and taking into account the equation ( pQ]) we get 
the relation 

^,|_(T,,./.)^o, ,3a) 

In our further considerations we neglect the interactions between the particles and 
suppose that the metric of the background gravitational field ga/s does not depend 
on motion of particles. After averaging ( p6D on paths we define the equation for the 
one-particle distribution function 

^-|;(W*)^o. ,3r) 

As a matter of principle we can generalize the problem |^| when particles are 
considered to be also sources of gravitational field which is self-consistent ly defined 
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from the Einstein equations with the energy-momentum tensor defined by the for- 
mula @. 

Taking into account the identity 

we get from (^) the colhsionless kinetic equation for the distribution function 

= 0, (38) 

where 

generahzes the Cartan's covariant derivation |^ for a space with higher order an- 
isotropy (locally parametrized by supporting elements (m",P/3)), provided with a 
(extended to momenta coordinates) higher order nonlinear connection 

Nl = 5l5lNt - T^P" (40) 

(on the geometry of higher order anisotropic spaces and superspaces and possible 
applications in physics see 0, [3l[]). 



The equation ( pH]) can be written in equivalent form for the distribution function 

Dxf = 0, (41) 

with the Cartan's operator defined on n(a) — 1 dimensional momentum space. 

The kinetic equations (|38|) and, equivalently, (^) refiect the conversation low of 
the quantity of particles in every volume of the space of supporting elements which 
holds in absence of collisions. 

Finally, in this subsection, we note that the kinetic equation ( |38D is formulated 
on the space of supporting elements {u°,Pj3) , which is characterized by coordinate 
transforms 

x'' = / {x') , / = y^K' {x') (42) 

and 

Ox^ 

Pj' = Pig— 7, Pa' = PaK^'ix'), 

where K^' (a;*) and K^,{x^) take values correspondingly in the set of matrices para- 
metrizing the group of linear transforms GL[m, JR"^), where IR denotes the set of 
real numbers. In a particular case, for dimensions n = m we can parametrize 

A distinguished by higher order nonlinear connection (^) tensor 
^Ikpl 'l^q i'^^iP-r) being contravariant of rang r and covariant of rang q satisfy the 
next transformation laws under chaingings of coordinates of type (^21) : 
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Operators of type (^) and transformations (^) and (53), on first order of anisotropy, 
on tangent bundle TM on a n dimensional manifold M were considered by E. Car- 
tan in his approach to Finsler geometry and by A. A. Vlasov |^6| in order to 
formulate the statistical kinetic theory on a Finsler geometry background. It should 
be emphasized that the Cartan- Vlasov approach is to be applied even in general 
relativity because the kinetic processes are to be examined in a phase space pro- 
vided with local coordinates {u'^^Pr)- Our recent generalizations 0, ^ to higher 
order anisotropy (including spinor and supersymmetric spaces) are to applied in the 
case of models with nontrivial reductions (modelled by N-connections) from higher 
dimensions to lower dimensional ones. 



4 Kinetic equation with pair collisions 

In this section we shall prove the relativistic kinetic equations for one-particle distri- 
bution function fiu^^pp) with pair collisions We summarize the related results 



and generalize the constructions for Minkowski spaces |T5 . 



4.1 Integral of collisions, differential cross— section and ve- 
locity of transitions 

If collisions of particles are taken into account (for simplicity we shall consider only 
pair collisions), the quantity of particles from a volume in the space of supporting 
elements is not constant. We have to introduce a source of particles C {x,p) , for 
instance, in the kinetic equations (^l]) , 

p"DJ = Cif) = Ciu,p). (44) 

The scalar function C {u,p) is called the integral of collisions (in the space of sup- 
porting elements). The value 

A"Wm —Lc(u,p) 

is the changing of the quantity of particles under pair collisions in a region 
Let us denote by 

W {p,P[i]\p',p\i]) 



P^P[i]P'^P[\] 



(45) 



the probability of transition for two particles which before scattering have the mo- 
menta pr- and p-Q and after scattering the momenta p^ and p'_p- with respective in- 
accuracies A"('')~V' and A"('')~Vp The function (the so-called 
collision rate) is symmetric on arguments p,P[i] and p',p[^ and describes the veloc- 
ity of transitions with conservation of momenta of type 
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when the conditions (^) hold. 

The number of binary colhsions within a la-spacetime interval 5u around a point 
u between particles with initial momenta in the ranges 

{lf,lf + S^) and ("p'[i],"p'[i] +(5"p'[i]) 

and final momenta in the ranges 

(jf', Y + 6Y) and ("p^Jij , "p^Jij + Slf[^) 

is given by 

f {u,p) f (u,p[i]^W (p, p[i] ; p', p'l]) (5?p5?p,^j5v'^^Pjij'^" (46) 

were, for instance, we abbreviated 5<jp = ^'"''•'^''p/ {\f\g\p'^)- 

The collision integral of the Boltzmann equation in la-spacetime (|4^ ) is expressed 
in terms of the collision rate (^) 

C if) = J f iu,p) f (u,piiyj (47) 

The collision integral is related to the differential cross section (see below). 



4.2 The Cross Section in La— spacetime 

The number n(^bin) of binary collisions per unit time and unit volume when the initial 
momenta of the colliding particles lie in the ranges 

Cf,^ + 6^) and 

and the final momenta are in some interval <^ in the space of variables Jf' and "p^j^j 
is to be obtained from (^) by dividing on Su = 6t6x'^...6x^dy^...dy"^ and integrating 
with respect to the primed variables, 

n(bin) = f {u, p) f {u, X 5^p5?-p[„ ^ W (p, p[i] ; p', p'^]) ^v^V^j • (48) 

Let us introduce an auxiliary velocity 

v = F/ip^pl) (49) 
were the so-called MoUer fiux factor is defined by 



— A 

m 



1/2 



(50) 



It may be verified that the speed v reduces to the relative speed of particles in a frame 
in which one of the particles initially is at rest. We also consider the product of the 
number density of target particles with momenta in the range "p^[i] + 

(given by / (^u,p[i]j S^^'^^^'^pm) and the fiux of incoming particles with momenta in 
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the range [jf, + Sjf) (given by / {u,p) v, where v is the speed 

product can be written 

Ff 5?p^%ij. 

By definition (hke in usual the relativistic Boltzmann theory 
section a is the division of the number (|48|) to the number (|5lD 



|T5|, 



a 



1 /" 

- y^l^(p,P[i];p',P[i])5v5s;i]- 



This 

(51) 
the cross 

(52) 



The condition that collisions are local implies that the collision rate must contain 
n(a) delta-functions 

W {p,p[i]]p',p[^) = w [p,piiY,p',p[^^) X + -p' -pjy) . (53) 

Let 1R"(-)-^ be a (^(a) — l)-dimensional Euclidean space enabled with an or- 
thonormal basis (^ei, 62, Cn^^j-i) and denote by (^m^, u^, the Cartesian 
coordinates with respect to this basis. For calculations on scattering of particles it 
is useful to apply the system of spherical coordinates (r„, , = r, 9n, ,-1, On, .-2, Oi) 

^ (a) (a) [a) / 

associated with the Cartesian coordinates {ui, U2, ■■■,) ■ The volume element can be 
expressed as 

where the element of solid angle, the (n(Q) — l)-spherical element, is given by 



dn+m-2 Q 



sm 



n+m—3 



On 



-2 ■ sm 



n+m— 4 



n+m—3 



sin^ ^3 ■ sin O2 x dOn 



+m-2 



.. ■ d02 ■ dOi. 



(54) 



In our further considerations we shall consider that the Cartesian and spherical 
(^(a) — 1) -dimensional coordinates are given with respect to a la-frame of type (^. 

In order to eliminate the delta functions from (^) put into (0) we fix as a 
reference frame the center of mass frame for the collision between two particles with 
initial momenta Jf and "^[ij- The quantities defined with respect to the center of 
mass frame will be enabled with the sub index CM. One denotes by 'j^cm ^^e polar 
axis and characterizes the directions 'p>^'cm of ^^e outgoing particles with respect to 
polar axis by means of generalized spherical coordinates. In this case 



P CM 



d 



P CM 



dVL 



CM 



(55) 



where dVLcM is given by the formula (^41) . 

The total {n{a) + -momenta before and after collision in la-spacetime are given 



by d-vectors 
Following from 



pc 



p P =P +P[i] 



(56) 



^CM 



'^PCM 



P CM 
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we have 



P CM 



P CM 



I ^CM '^^CM- 



In result the formula for the volume element fpSl) transforms into 



P CM 



1(a)- 



P'cwi dPcM dVtcM- 



Inserting this collision rate ( p3D into (|5^ ) we obtain an integral which can be 
rewritten (see the Appendices 1 and 7 to the paper Q for details on transition to 
the center of mass variables; in la-spacetimes one holds good similar considerations 
with that difference that me must work with respect to la-frames (^ and (|iy) and 
d-metric (0) ) 



CM 



(57) 



where E = y/P^Pa, , is the total energy (devided by c = 1). 

The differential section in the center of mass system is found from ( p?] ) 



6a_ 



I CM 



which allow to express the collision rate 

W {p,P[i];p',p[i] 
We shall use the formula 



_P"'(a)-4 

as 



w 



I CM 



P[l] -p - P[l] 



in order to compute the transport coefficients. 



(5^ 



5 Equilibrium States in La— Spacetimes 

When the system is in equilibrium we can derive an expression for the particle 
distribution function f {u,p) = f(^eq) {u,p) in a similar way as for locally isotropic 
spaces and write 

/m P) = ^2nnn^-^ [ j ' 

where is Boltzmann's constant and h is Planck's constant divided by 27r and = 
fi {u) and T = T {u) are respectively the thermodynamic potential and temperature. 
For simplicity, in this section we shall omit explicit dependencies of fi and T on la- 
spacetime coordinates u. Our thermodynamic systems will be considered in local 
equilibrium in a vicinity of a point U[o] with respect to a rest frame locally adapted to 
N-connection structure, like (|^) and (p^Of). In the simplest case the n + m splitting is 



trivially given by a N-connection with vanishing curvature (|TT|). For such conditions 
of trivial la-spacetime and T can be considered as constant values. 
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5.1 Particle Density 

The formula relating the particle density n, temperature T, and thermodynamic 
potential /x is obtained by inserting (|59D into (pTD , with (^71) . Choosing the locally 
adapted to the N-connection to be the rest frame, when = (1, 0, . . . , 0) the 
calculus is to be performed as for isotropic n(a) dimensional spaces The integral 
for the number of particles per unit volume in the rest la-frame is 

n{,^^T) = \ .exp ( X / /exp (-^] (60) 

with 

/= (|p^|2 + ^2^^/^ (61) 

and depends only on length of n(a) — 1 dimensional d-vector "p^. By applying spherical 
coordinates, when 

and dVt is given by the expression (p4|), and differentiating on radius p the well 
known formula for the volume Vn^^^^-i in a n(^a) ~ 1 dimensional Euclidean space 



(rz (,) - l) vr("(")-0 /^ 
r[(n(,) + l)/2] 



^-w-i = HtfTI-^^YT^A (62) 



where F is the Euler gamma function, then putting p = 1 we get 

/ = ^ Zu ^ ^T^- (63) 

J r[(n(,) + l)/2] ^ > 

We note that from (|6TD one follows 

and we can transform (^) into an integral with respect to the angular variables 
through the replacement 

pT r|(n,„, + 1)/21 

Introducing the dimensionless quantities 

^ = p^/ksT, and ^9 = m/keT, (64) 

for which, respectively. 



rfpi = A;bT di, and | p^l = A;BrA/^2 _ (65) 



22 



the integral (^) is computed 



ni,,T) ^ 2"-/^ ^ (66) 



where the modified Bessel function of the second kind of order n(a)/2 has the integral 
representation 



We note the dependence of (|66|) (but not of (|65D ) on m anisotropic parame- 
ters (coordinates). The formulas proved in this subsection transforms into locally 
anisotropic ones if the N-connection is fixed to be trivial (with vanishing N- 
connection curvature (|11|)) and the d-metric (|I3|) transforms into a usual (pseudo) 
Riemannian one. 



5.2 Average energy and pressure 



The average energy per particle, for a system in equilibrium, can be calculated by 
introducing the distribution function (|59|) into (|33|) and applying the formulas (pSf), 
and (p^). In terms of dimensionless variables (B^ we have 



{kBTf^''^ (nia) - l) 7r("w-^)/' 



(27r;i)"("' 
/i 



r[(^{a) + 1)/2] 



X exp 



kfiT J J'& 



("(a)-3)/2 



After carrying out partial integrations together with applications of the formula 



(s-2)/2 _l d 

s dq 



s/2 



we obtain 



TT 



X exp 



(27r;i)"w-' UJ 

^ ^ ^^^(n(.)+2)/2 (^) - i^n(„,/2 (^); 



(68) 



Dividing the energy density (|68|) to the particle density (|60D and substituting 
the inverse dimensionless temperature i} = i} (T) (^) we get the energy per particle 



e (T) = ksT 



m ^{n(„)+2)/2 
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(69) 



This formula is the thermal equation of state of a equilibrium system having m 
anisotropy parameters (we remember that n(^a) = n + m) with respect to a la-frame 

®- 

The pressure d-tensor is to be computed by substituting (^) into (p4D, applying 
also the integral (pSl) . We get 

p^P = -keT n in, T) A"^, (70) 

were, by definition, the coefficients before determine the pressure 

P = kBT nifi,T). (71) 

So, for a system of n(^a) — l=n + 'm — 1 dimensions the expression ( [7TD defines 
the equation of state of ideal gas of particles with respect to a la-frame, having m 
anisotropic parameters. 

5.3 Enthalpy, specific heats and entropy 

The average enthalpy per particle 

h{T) = e{T)+P{fi,T)/h{fi,T) (72) 

is computed directly by substituting in this formula the values (H),(0) and (|66D . 
The result is _ 

t(r)^a ^^y-';-(ff . (73) 

Using (^) and (^) we can compute respectively the specific heats at constant 
pressure and at constant volume 

Subtracting of h and e, after carrying out the differentiation with respect to T, 
one finds Mayer's relation, 

and, introducing the adiabatic constant 7 = Cp/c„, 



The relation (|7^ can be proven by straightforward calculations by using the prop- 
erties of Bessel's function (|67|) (see a similar proof for isotropic spacetimes in [§). 
The entropy per particle is introduced as in the isotropic case 

. = (75) 
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We have to insert the values of fi (found from 

/i (n, T) = kBT X In (^Tr^fw |n 

(2m)«(^)/2(vrfc^T)K^)-2)/2ir„^„,/,(^) 

and of h (see (|73[) in order to obtain the dependence s = s (T) (for simphcity, we 
omit this cumbersome formula). 

5.4 Low and high energy Umits 

Let us investigate the non-relativistic locally anisotropic limit by applying the 
asymptotic formula [|1| 

valid for large ^. The formula (^) gives 

~ n+m-1 
e m-\ kBi . 

Thus, each isotropic and anisotropic dimension contributes with ksT /2 to the aver- 
age energy per particle. 

Now, we consider high temperatures. For t9 — > one holds the asymptotic 
formula [0] 

Ks {^) ~ 2'-^T (s) 

and the particle number density (|5BD, the energy (|5U[), enthalpy ( [73| ) and entropy 
(|75D can be approximated respectively (we state explicitly the dimensions n + m) 

by 

/ I \ / h rp\n+m-l 
-,n+m-ln+m-2-r / ^ + \ / I^BJ- \ ( 



fi T) = 2— V»--r j X j exp J . (77) 

e = {n + m — 1) ksT, h = [n + m) ksT, 
s = {n + m)kB- 

In consequence, the corresponding specific heats and adiabatic constant are 

n + m 



{n + m) ks, Cy = {n + m — 1) ks and 7 



n + m — 1 



These formulas imply that for large T the spacetime anisotropy (very possible at 
the beginning of our Universe) could modify substantially the thermodynamic pa- 
rameters. 

Finally we emphasize that in the ultrarelativistic limit 



,.r-..--,„.,„,xr(!i±^)(|I)"^'"-'|. ,78) 

We also note that our treatment was based on Maxwell-Boltzmann instead of Bose- 
Einstein statistics. 



25 



6 Linearized Locally Anisotropic 
Transport Theory 

For systems with local anisotropy outside equilibrium the distribution function 
f{eq){u,p) (|59|) must be generalized to another one, f{u,p), solving the kinetic la- 
equation (PI). We follow a standard procedure of linearization and construction of 
solutions of kinetic equations by generalizing to la-spacetimes the Chapman-Enskog 
approach (we shall extend to locally anisotropic backgrounds the results presented 



m 



We write 

fiu,p) = f[o]iu,p)[l + ifiu,p)] (79) 



with the lowest order of approximation to / taken similarly to (59) 



where the constant variables of a trivial equilibrium fi, T and Ua are changed by 
some local conterparts fi{u) ,T (u) and Ua (u) . Outside equilibrium the first two 



dependencies /i (u) and T (u) are defined respectively from relations ([Tq) and (69). 



6.1 Linearized Transport Equations 

The so-called deviation function ip {u, p) from (|79D describes the deformation by 
nonequilibrium flows of /[oj (m,p) into / {u,p) . It is considered that in equilibrium ip 
vanishes and not too far from equilibrium states it must be small. The Chapman- 
Enskog method states that after substituting /p] into the left side and /pj (1 + v?) 



into the right hand of (Q) we shall neglect the quadratic and higher terms in (p. 
In result we obtain a linearized equation for the deviation function (for simplicity, 
hereafter we shall not point to the explicit dependence of functions, kinetic and 
thermodynamic values on spacetime coordinates) 

-p"^a/[0] = f[Q\LM 

where, having introduced (|80D into (0), we write for the right side 

' H-p'Ui -p'^Ua\ 



The left side is to be computed by applying the generalized Cartan derivative ([391). 

We introduce a locally anisotropic generalization of the gradient operator by 
considering the operator (in brief, la-gradient) 

V„ = Dp (81) 

with given by (^) (we use the a d-covariant derivative defined by a d-connection 
([mD instead of partial and / or isotropic derivatives for isotropic spaces [l^, |^ ) . In a 
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local rest frame of the isotropic fluid Va (0, di, 82, ■■■) , i.e. in the flat isotropic 
space the la-gradient reduces to the ordinary space gradient. By applying Vo we 
can eliminate the time la-derivatives 

knTLl^] = py^V^U, - Tp-V^ (^) - pVf/« - ^) (82) 



+[(7 - i)p''u^ + ^ - 1) It j + ^ 1,9^ jVu.v^u^. 

We can verify, using the expression /i = /i (n, T) from (^) that one holds the 
equalities 

d / fi\ e (9/i ksT 
df \f) ~ ^r2' d?i ~ 

and 

Introducing these expressions into (p2D we get 



where 



ksT L[^] =EX + {h- p^U^) A^^ppX" 

+ (a^ A^; - -—^ A,,A^") p0^ 

V ' n + m — 1 



, -7 (p°f/a)' + [(7-l)/^-7A:BT]p"f/„-— (83) 

.n + m — 1/ n + m — 1 

and there were considered forces deriving the system towards equilibrium 

V" T V° » 
X = -V'U., X" = -—- - --^ 
^ T nh 

and 

(o)X"^ = - N'W + V^t/n ^ A^^V^f/.. 

2 ^ ^ n + m - 1 

We note that the values ksT L[ip\, S and ^^'^X^^ depend explicitly on the dimen- 
sions of the base subspace, n, and of the fiber subspace, m. The anisotropy also 
modify both the thermodynamic and kinetic values via operators and A^ which 
depend explicitly on d-metric and d-connection coefficients. 

6.2 On the solution of locally anisotropic transport equa- 
tions 

Let us suppose that is known the solution of these three equations: 

kBTL[A{p)] = E, (84) 
kBTL[B{p)A^^Pp] = {h-p-U^)Aipp 
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and 



keTL 



C (P) (A^A,- - 



1 



n + m — 1 
1 



A„,A^'^^ p^p^ 



n + m — 1 

for some functions A{p) ,B {p) and C {p) . Because the integral operator L is linear 
the linear combination of thermodynamic forces 



ip = AX + + Cap (°)X"^ 



where 



and 



Ba = Bip) /\i pp 



Cap = B{p)( Al^AZ - 



-AapA' 



Pe Pa, 



n + m — 1 

is a solution for the deviation function ip. 

It can also be verified that every function of the form if = a + baP"' with some 
parameters a and ba not depending on is a solution of the homogeneous equation 
L [ip] = 0. In consequence, it was proved (see and 0) that the scalar parts of A (p) 
and B (p) are determined respectively up to functions of the forms (f = a + baP" 
and b°'pa- We emphasize that solutions of the transport equations (|8^ (see next 
subsections) will be expressed in terms of A, B, and C. 



6.3 Linear laws for locally anisotropic non— equilibrium 
thermodynamics 

Let Up{u) be a velocity field of Landau-Lifshitz type characterizing a locally an- 
isotropic fluid flow. The heat flow and the viscous pressure d-tensor with respect to 
a such la-field are correspondingly defined similarly to but in terms of objects 
on la-spacetime (see formulas (|T3|), (^, (|3|), (If, (^, (^, (|g), (0), (|TUP, (|g) 



liHeat) = [UaT-^ - hn^^) A^^, (85) 
Yjal3 _ paf3 _|_ p/\al3 



It should be noted that if the Landau-Lifshitz condition (30) is satisfied the first term 
in I^heat) vanishes and the heat flow is the enthalpy carried away by the particles. 
The pressure P (u) is defined by n{u) ksT {u) , see ([71]). Inserting (^) and (|80|) 
into (^) we prove (see locally isotropic cases in |T5|, ^) the linear laws of locally 
anisotropic non-equilibrium thermodynamics 

^{heat) = ^TX" 

and 
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with the transport coefficients (the heat conductivity A, the shear viscosity r], and 
the volume viscosity coefficient ?7(^); in order to compare with formulas from we 
shall introduce explicitly the light velocity constant c) defined as 

(n + m — 1)1 J iji- 



n + m — 1)T 

c 

n + m 
(n + m) (n + m — 4:) J 



V{v) = 7 / PePaA^" Afyo^dqp, (86) 

^ ' n + m — 1 J I I 



{n + m) {n + m — A) 

where one uses the d-tensor 

1 , . , 1 



A'i't = - (AlA^ - A^A^t) A'^^A, 



with the properties 
and 



A pa ^ _i , in + m){n + m-l) 

^pa 2 



The main purpose of non-equilibrium thermodynamics is the calculation of trans- 
port coefficients. With respect to la-frames the formulas are quite similar with those 
for isotropic spaces with that difference that we have to consider the values as d- 
tensors and take into account the number m of anisotropic variables. 

On both type of locally isotropic and anisotropic spacetimes one holds the so- 
called conditions of ffi (see, for instance, 0) 

/ PeU'Af[o]d^p = and J {p.U'f Af^.^d^^.^ = 

which allow us to write the volume viscosity (see (p3D) 

?7(„) = c j EAf[o]dqp. 

For some d-tensors i;/""!'---"? (^^^ a^d Sai,...aq{p) "we define the symmetric bracket 

{H, S} = ^J if"--"' (p) L [S^,,...aM] fioAp 

where n is the particle density and L is a linearized operator. In terms of such 
brackets the coefficients (|86|) can be rewritten in an equivalent form 

A = - (87) 
n + m — 1 



(n + m — 1)^ — 3 
r]{v) = cksTn^ {A, A} . 

The addition of invariants of type ip = a + b^p"' to some solutions for A (p) 
Ba (p) , or Caf3 (p) does not change the values of the transport coefficients. 
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6.4 Integral and algebraic equations 

The solutions of transport equations ( p^ are approximated by considering power 
series on ( = p^Ua/ksT for functions 

oo oo ^ oo 

^ (C) = E ^aC, ^ (0 = E ^ (0 = E c^^- 

a=2 "6=1 c=0 

The starting values 2 = 2 and 6 = 1 have been introduced with the aim to determine 
the scalar functions A {Q and B {Q , respectively, up to contributions of the forms 
a + 6oP" and h'^Pa- Inserting these power series into the integral equations (|5^, 
multiplying respectively on C*/[o]! CP^/p], and Cp'^p'^fio], after integrating on d^p 
we get 

oo oo oo 

E «ai a^a = ' E fe^ = ' ^ci^C*- = 7?^ , 
a=2 b=i c=0 

where Si = 2, 3, . . . , 61 = 1, 2, . . . , ci = 0, 1, . . . , and there are symmetric brackets 

= {C\C}, (88) 

= {CV,^AX}, 



c- 



cic 



and integrals 



C- (ft - p °I/„) A ^.p^pVioi 



, ^TpPuPf,p"p^ f[0] ^ 



The lowest approximation is given by the coefficients 

^2 = "2/ «22' % = /^/^' = V Coo- 
Introducing these values into (^) we obtain the first-order approximations to the 
transport coefficients 

A = (90) 

n + m — 1 



+ m) (n + m — 1) — 2 c^' 



?7(„) = cksTn^^. 

^22 
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The values a^, /J^^ and are (n + m — l)-fold integrals which can be expressed in 
terms of enthalpy h and temperature T, 



7n 



n + m + 1 — (n + m — 1)7 (n + m — 1)7 
/c^T (n + m — 1)7 ksT 



0, 



cn 



7-1 



cn 



cn 

T^oo 

(n + m — 1) (n + m) 
2 



(91) 



^ra + ni — 2) (n + m + 1) /i 



1 kuT 



n \ c 



in + m) — 1 



(n + m) 



In order to complete the calculus of the transport coefficients it is necessary to 
compute the brackets from the denominators of (pO|) . 



6.5 Brackets from flat to la— spacetimes 

The calculation of brackets is a quite tedious task (see |^) which should involve the 
curved spacetime metric and connection. For simplicity, we consider a background 
flat spacetime with trivial local anisotropy. In this case we can apply directly the 
formulas proved for the Boltzmann theory in n{^a) dimension but with respect to 
la-frames and by introducing the locally anisotropic d-connections and d-metrics 
instead of their isotropic analogous. 

Let us consider a locally adapted to N-connection rest frame, ~^ being a unit 
vector and denote by (^9„j^j„i, 6„^^^_2, 61^ the angles of spherical coordinates 
with respect to a Cartesian one (X^,X^, chosen as the X'^('")-axis to be 



along 'fcM (where \'fcM\ = F/P,P = \P" 



see 



Oj) (ill), and if CM ~^ ™^ ^^e 



plane of X^C") and axes. We suppose that the scattering angle is 6„j^j_i and 

the differential cross section {da/dfl)(-,j^^ is a function only on the module of total 
energy and on scattering angle 6 = 6„^^j_i, 



da' 

dn 



CM 



' da 

dn 



CM 



The results of calculations of symmetric brackets will be expressed into terms of 
twofold integrals, which for arbitrary integers 



J- 



X 



1 



r 



s \ ( d/2\ 2^+""'^ 



t j \ u 

00 



TT 



(92) 



F(t,u,v) / dx X^" 



i+m—l+d—u—b—v+i) 



2C 



/ 2 \ {n+m-2+t+u+j)/2 



■'Ow 



U 
W 



sin"' e cos"-'" e 



X 



2C 



CM 
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where there are introduced the integers 

d = r + s-t-u,b = (t + u + n + m-2^ /2,/i= [1 - (-1)'^] 
and variables x = cj (ksT) and ( = mc^ / (ksT) . The factors in (p^ ) are defined 

^(?i+m-3)/2 



HO 



2T ((n + m - l)/2) T (n + m - 3) C Kn+rn {Q 



2 1 



and 



F (J, id^ 



l + i-l) 



l + i-l 



.t+u—w 



, n + m — 3 id + 1\ ( n + m — 2 + w t + u — w + 1 
xB , — - — \ X B { — 



2 ' 2 y V 2 

where the beta function is given by gamma functions 

r(x)r(y) 



B {x,y) 



and by f - j it is denoted the Newton's binomium. 



6.5.1 Scalar type brackets 

The first type of brackets necessary for calculation of transport coefficients (|87|) (see 
the series approximation ( pH] ) are the so-called the scalar brackets, decomposed into 
a sum of two integrals 

«aia = <a + ^^a' (93) 

where Si, a = 2, 3, 



J, . . . , 



aia 



n ^ j f[o] {p)f[o] {p[i]) iCpT' X {{Cpi^y - (Cp;J1w^f^?"pf^^P[i]C?V^%] 

^~^/ /[o] ip)fm {p[i]) iCpT' X liCpT - (Cp')"]w^f^s^^P[i]*P'^s;i]' 

with the integration variables 

Cp = p''UjkBTXpy,,=P[ifJkBT, 
Cp' = p"'UjkBTXp'=p'(^]UjkBT. 



The explicit calculations of integrals from (|93|) (see the method and basic intermediar 
formulas in P]; in this paper we deal with la-values) give 

ai a d/2 u 



a\ a u 
4=0 «=0 v=0 w=0 
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In the lowest approximation, for Si = a = 2, we obtain 



'^22 ~ ^ 



J|'°'°) (2,2,0,0)+j|°'°) (2,2,0,2 



So, we have reduced the scalar brackets to twofold integrals which are expressed in 
terms o spherical coordinates with respect to a locally anisotropic rest frame. 

6.5.2 Vector type brackets 

The vector type brackets from (|87|) and (|8^) are also split into two types of integrals 

^.1^=% + % (94) 

where 6i, 6 = 1, 2, . . . , 



X 



Wdqpdqp.^Aqp'd^p 



bib 



A tedious calculus similar to that presented in implies further decompositions of 
coefficients and their representation as 

3 3 

bib ^ {t)bib bib ^ (r)bib^ 



(i)=l 



with corresponding sums 



b' 



(i)feife 



4 \ c 



X 



EEEE(-ir4r {I 



U, V, w 



h' 



t=0 11=0 v=0 10=0 
2 bi 6 d/2 M 



(2)bib 



b', 



(2)bib 



i-f) 1:1:1:1: 

^ / t=0 u=0 v=0 w=0 

/ U rpX'i bi b d/2 u ^ 

- m i:i:i:i:(-ir4;ii,( 

^ / 4=0 u=0 v=0 10=0 



with corresponding sum expressions for b"^.^^ by omitting the factor (—1)" . 

In the lowest approximation (for 61 = 6 = 1) we have 

2 



b^ 



X 



[4?'"'^ (1,1,0,0 



+J|'°'°) (2,2,0,0)+j|'°^°) (2,2,0,2 



Here should be noted that in general T = T (u) is a function on la-spacetime 
coordinates. 
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6.5.3 Tensor type brackets 

In a similar fashion as for scalar and vector type symmetric brackets from (^) and 
(|88D one holds the decomposition 

Ccic = 4c + c^ic (95) 

where ci, c = 1, 2, ... , 

c^T^c = / fm ip) /[o] {pm) K'e iCpf p.p. 



ClC 



x[Pli]P[i] (Cp[i,)'-p'T]P'i] (Cp|^j)']Vr*prf?pjyC/v*) 

= n-' f fmip)fm{pii])A:^ACpfpuP, 



X 



"11] J J '-f -'-lij 't- 'I'll] 
For tensor like brackets we have to consider sums on nine terms 

c^-= y c;,--andc^-= y c''--. 

ClC / ^ (•5)ClC ClC / ^ (S)CIC 

(s)=l (s)=l 

These terms are (with that exception that we have dependencies on the number 
of anisotropic variables and T = T (m) is a function on la-spacetime coordinates) 



^ / t=Q u=0 v=0 w=0 



'knT 



4 ci c d/2 u 



= Ff EEEE(-ir4r('. 

^ ^ t=0 «=0 ■u=0 »i)=0 

./ _ n + m-2 (ksTV ^J^^ ^ ^ ^ « .(o,o,o) 



'(3)cic ^ _^ ^ _ X I c 



EEEE(-i)"4S2(f 



C.oNCCr = I I > > > > (-i) , 



t=0 M=0 D=0 10=0 



1 /fcpT 



4 ci c u 



= EEEE(-ir4r'('. 



t=0 «=0 j;=0 ui=0 



4 ci c d/2 M 
C' - - - O / ^ A" 7(0.2-1) 



EE EE(-ir4;,L 



{5)cic \ C I ^ ^ ^ ^ ci+lc+1 

^ ' ^=0 «=0 ?5'=0 10=0 



^ ^ t=0 Jl=0 ?;=0 10=0 



/'4 ( h. 4 ci c rf/2 M 

- -^TT^pf SEE i:(-i)"4r 

^ / t=0 n=0 D=0 i«=0 



A4 / 1, \ 4 ci ^ (i/2 n 

v.. = ^rri^ m EEEE(-ir4S(?.fi.--.'"). 

^ ' t=0 u=0 v=0 10=0 
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The twice primed values c"^-,^^ are given by similar sums by omitting the factor 
(— 1)" and by changing into c"^j^^-^ and c'^^j^^-^ the overall signs. 



In the lower approximation, for Ci = c = 0, one holds 

4 



knT' 



n + m 



{Ji^'^'^^(0,0,0,0) + 2Jii!'^'^^(l, 1,0,0) 



11 



,„ + ,„_ll4''°*P.2.5.0) + jt°-»'(2.2.0.2)]}. 

6.6 Locally anisotropic transport coefficients for a class of 
cross sections 



For astrophysical applications we can substitute 



da 



(96) 



CM 



with some scalar factor ^ (with or without dimension) and r being a positive or 
negative number into (|92D . In this subsection we put ( = mc^/ [ksT) for high 
values of T. The chosen type of differential cross section (p6| ) is used, for instance, 
for calculations of neutrino-neutrino scattering (when r = 2 and ^ is connected with 
the weak coupling constant). In the first approximation the symmetric brackets ( P5| ) 
are 



"22 



hi 



"00 



'2kBT\ 7r('^+™-i)/2r [n + m + l + r/2]T[{n + m + r) /2] 
c J T[n + m-2]T[{n + m) /2] T[{n + m+l) /2] ' 

/ kBT\^ n + m + 2 + r 



(97) 



'knT' 



2 22' 

(n + m + r) (n + m + 4 + r) 



+ 



n + m 



n + m — 1 



"22- 



Putting these values into (|9T1) we obtain the locally anisotropic variant of trans- 
port coefficients in the first approximation, when ?7(^) ~ but with nonzero 



A 



2/cbc / c \ (n + m)^ (n + m — 1) 
\2kgr) n + m + 2 + r 
T [n + m - 2]T [{n + m) /2] T [{n + m + 1) /2] 
T[n + m + l + r/2]T[{n + m + r) /2] 



(98) 



and 



e \2kBT 



r-l 



^{n+m-l)/2 



(n + m) —12 — 171 — 2 
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^ [{n + m-lf-l\ (n + m + 1) ^^^^ 

(n + m + r) (n + m + 4 + r) + 2 ["^^I^j 

r [n + m - 2] r [(n + m) /2] T [(n + m + 1) /2] 
^ r + m + 1 + r/2] r [(ra + m + r) /2] 

These formulas present a locally anisotropic generalization of the Boisseau and 
van Leeuwen [01 results. 



7 Transport Theory in Curved Spaces with Rota- 
tion EUipsoidal Horizons 

After having established in a general way the scheme for calculation of transport 
coefficients on la-spacetimes, we now specify an example for a four dimensional 
static metric (being a solution of Einstein equations in both general relativity and 
la-gravity) with the event horizon described by a hypersurface of rotation ellipsoid 
(see the Appendix). We note that in this case n = 3, m = 1 and n<^a) = 4. 
The main formulas for kinetic and thermodynamic variables from the previous 
Sections 3-6 were based on spherical symmetry of n(a) — 1 dimensional volume (|62|) 
and spherical integral (§3|). For the rotation ellipsoid we have to modify the volume's 
formula by introducing the ellipsoidal dependence 







1 — 

r 


\n(a) - l)/2 





Hay 



=4 



were 

^ ^ 1 -ecos^ 

is the parametric formula of an ellipse with constant parameter p(o), angle variable 
6 and eccentricity e = 1/a < 1 is defined by the axes of rotation ellipsoid (see the 
formula ( |119D from Appendix). 

If in the case of spherical symmetry 4 J d6 = 27r, the ellipse deformation gives 



the result 

tt/2 



, , dO 8 l-e 
4 / ^ = ^^^^= arctan W . 







So, performing integrations on solid angles in spaces with rotational ellipsoid sym- 
metry we can use the same formulas as for spherical symmetry but multiplied on 



4 l-e 
= — arctanW^— . 



For instance, the integral (B^) transforms 



(n(.) - l) 7r(»(-)-^)/^ 3vr3/2 
Ha) - l)/2 i [<^/2J 
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The formulas for the particle density ( |5DD and energy density (^) of point par- 
ticles must be multiplied on 

n n(e) = q(e)n and e = q(e)e, 

but the energy per particle will remain constant. We also have to modify the formula 
for pressure (|7TD, been proportional to the particle density, but consider unchanged 
the averaged enthalpy ([7^). The entropy per particle ([75|) and chemical potential 
( |76D depends explicitly on g(j)-factor because theirs formulas were derived by using 
the particle density n(^^y Here we note that all proved formulas depends on m(= 1, 
in this Section) anisotropic parameters and on volume element determined by d- 
metric. In the first approximation of transport coefficients we could chose a locally 
isotropic background but introducing the factor g(£) and taking into account the 
dependence on anisotropic dimension. 

Putting e-corrections into (pO| ) we obtain the first-order approximations to the 
transport coefficients in a la-spacetime with the symmetry of rotation ellipsoid 

A = (100) 







10 c^- 

^22 



The values a^, and 7^ are 3-fold integrals expressed in terms of enthalpy h and 
temperature T and have the limits 

«2 = ^C^O' (101) 



To 
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ch^e) 7 - 




In the next step we compute the -deformations of brackets from the de- 
nominators of (|100| ). Because the squares of a, (3 and 7 coefficients from (|89D are 
proportional to deformations of fi~^ (this conclusion follows from the formulas (|97|) 
in the T — >• cxo limit, see ( |101| )) and the scalar (0), vector (p^ ) and tensor (p5| ) 
type brackets do not change under g(e)-deformations (see (^)) we conclude that 
the transport coefficients ( pUD do not contain the factor but depends only on 
the number m of anisotropic dimensions. This conclusion is true only in the first 
approximation and for locally isotropic backgrounds. In consequence, the final for- 
mulas for the transport coefficient, see (^) and (p9D , in a (3+1) locally anisotropic 
spacetime with rotation ellipsoid symmetry are 
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2kBC / c y 48 r[2]^r[5/2] 

\2kBT) 6 + r r [5 + r/2] T [(4 + r) /2] ' 

1 / c y-i 20 r[2]^r[5/2] 

i \2k^) 7r3/2[(4 + r)(8 + r) + 4/3] T [5 + r/2] T [(4 + r) /2] ' 

This is consistent with the fact that we chose as an example a static metric with a 
local anisotropy that does not cause drastic changings in the structure of transport 
coefficients. Nevertheless, there are g(£)-deformations of such values as the density 
of particles, kinetic potential and entropy which reflects modifications of kinetic and 
thermodynamic processes even by static spacetime local anisotropies. 



A = 

V = 



8 Concluding Remarks 

The formulation of Einstein's theory of relativity with respect to anholonomic frames 
raises a number of questions concerning locally anisotropic field interactions and 
kinetic and thermodynamic effects. 

We argue that spacetime local anisotropy (la) can be modeled by applying 
the Cartan's moving frame method with associated nonlinear connection (N- 
connection) structures. A remarkable fact is that this approach allows a unified 
treatment of various type of theories with generic local anisotropy like generalized 
Finsler like gravities, of standard Kaluza-Klein models with nontrivial compactifica- 
tions (modelled by N-connection structures), of standard general relativity with an- 
holonomic frames and even of low dimensional models with distinguished anisotropic 
parameters. We have shown a relationship between a subclass of Finsler like metrics 
with (pseudo) Riemannian ones being solutions of Einstein equations. 

This paper has provided a generalization of relativistic kinetics and nonequilib- 
rium thermodynamics in order to be included possible spacetime local anisotropies. 
It should not be considered as a work on definitions of some sophisticated theories 
on Finsler like spaces but developing an approach to kinetics and thermodynamics 
in (pseudo) Riemannian spacetimes of arbitrary dimension provided with anholo- 
nomic frame structures. The item of introducing into consideration of generalized 
Finsler and Kaluza-Klein spaces was imposed by the facts that one has recently 
constructed (see Appendix and solutions of the Einstein equations with generic 



anisotropy, of Finsler and another type, like ellipsoidal static black holes, black tora, 
anisotropic solitonic backgrounds and so on) and that in the low energy limits of 
string theories various classes of generalized Finsler-Kaluza-Klein metrics could be 
obtained alternatively to the well known (pseudo) Riemannian ones By ap- 



plying the moving frame method it is possible to elaborate a general schema for 
defining of physical values, basic equations and approximated calculations of kinetic 
and thermodynamic values with respect to la-frames in all type of the mentioned 
theories. 

The crucial ingredient in definition of collisionless relativistic locally anisotropic 
kinetic equation was the extension of the moving frame method to the space of 
supporting elements (m",p'^) provided with an induced higher order anisotropic 
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structure. The former Cartan-Vlasov approach proposing a variant of sta- 

tistical kinetic theory on curved phase spaces provided with Finsler hke and Car- 
tan N-connection structures, was self-consistently modified for both type of locally 
isotropic (the Einstein theory) and anisotropic (generalized Finsler-Kaluza-Klein) 
spacetimes with N-connection structures induced by local anholonomic frame or via 
reductions from higher dimensions in (super) string or (super) gravity theories. The 
physics of pair collisions in la-spacetimes was examined by introducing on the space 
of supporting elements of (correspondingly adapted to the N-connection structure) 
integral of collisions, differential cross-sections and velocity of transitions. 

Despite all the complexities of definition of equilibrium states with generic anisot- 
ropy it is possible a rigorous definition of local equilibrium particle distribution func- 
tions by fixing some anholonomic frames of reference adapted to the N-connection 
structure. The basic kinetic and thermodynamic values such as particle density, 
average energy and pressure, enthalpy, specific heats and entropy are derived via in- 
tegrations on volume elements determined by metric components given with respect 
to la-bases. In the low and high energy limits the formulas reflect explicit depen- 
dencies on the number of anisotropic dimensions as well on anisotropic deformations 
of spacetime metric and linear connection. 

One can linearize the transport equations and prove the linear laws for locally 
anisotropic non-equilibrium thermodynamics. It has also established a general 
scheme for calculation of transport coefficients (the heat conductivity, the shear 
viscosity and the volume viscosity) in la-spacetimes. An explicit computation of 
such values was performed for a metric with rotation ellipsoidal event horizon (an 
example of spacetime with static local anisotropy), recently found as a new solution 
of the Einstein equations. 

Our overall conclusion is that in order to obtain a self-consistent formulation 
of the locally anisotropic kinetic and thermodynamic theory in curved spacetimes 
and calculation of basic physical values we must consider moving frames with cor- 
respondingly adapted nonlinear connection structures. 
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Appendix: A Locally Anisotropic Solution of Ein- 
stein's Equations 

Before presenting an explicit construction of a four dimensional solution with 
local anisotropy of the Einstein equations (^) we briefly review the properties of four 
dimensional metrics which transforms into (2+2) or (3+1) anisotropic d-metrics (we 
note that this is not a (space + time) but a (isotropic+anisotropic) decomposition of 
coordinates) by transitions to correspondingly defined anholonomic bases of tetrads 
(vierbeins). 
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Let us consider a four dimensional (in brief, 4D) spacetime (with two 

isotropic plus two anisotropic local coordinates) provided with a metric of signature 
(+,+,-,+) parametrized by a symmetric matrix of type 



gi + qi'^hs + ^1^/14 qihs nihi 

92 + g2^/i3 + n2^h4, q2h3 712/14 

qihs q2h3 

nihi n2hi /14 



(102) 



with coefficients being some functions 

Qi = gi{x^),qi = qi{x^,t),ni = ni{x^,t),ha = ha{x^,t) 

of necessary smooth class. With respect to a la-basis ( p!0D the ansatz (|102|) results 
in diagonal 2 x 2 h- and v-metrics for a d-metric (|T3|) (for simplicity, we shall 
consider only diagonal 2D nondegenerated metrics because for such dimensions every 
symmetric matrix can be diagonalized) . 

An equivalent diagonal d-metric of type (|13]) is obtained for the associated N- 
connection with the coefficients being functions on three coordinates {x'^,z), 



qi{x\z), = q2ix\z), 
ni{x\z), N2 = n2{x\z). 



(103) 



For simplicity, we shall use brief denotations of partial derivatives, like a = da/dx^, 
a' = da/dx'^, a* = da/dz a! = d'^a/dx^dx'^, a** = d'^a/dzdz. 

The non-trivial components of the Ricci d-tensor (^) (for the ansatz |102| ) when 
Rl = R2 and = Sf, are computed 



Rl 



Si 



R 



31 
P32 

P42 



Rr, 



q4 



^1 

2 ' 

-I 
2 ' 



1 . 



1 
2^ 



92 + 9i92 



1 

2^ 



/ 2 

9i 



9i92 



2 



h 



-K* 



h 

1 



2h? 



h*h* 
"'■i"'A 



2/13/14 2/l4 2/l4 



[^K -hi + 



2h. 



-h 



2/142 



2/13/14 



1 h' 



(104) 
(105) 

(106) 



h' 

h' * 4- — 



4 I 



(107) 



h*n — ?)h*,]nl 



The curvature scalar R (|T^ is defined by the sum of two non-trivial components 
R = 2R\ and S = 2S'|. 
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The system of Einstein equations 

R\ = 



transforms into 



SI 



KT4j, 



(108) 
(109) 
(110) 
(111) 



where the values of R\, S'f, Pai, are taken respectively from ( |104D , ( |105D , ( |106D , ( IIOT] ). 

We note that we can define the N-coefficients ( |103| ), qi{x^,z) and ni{x^,z), by 
solving the equations ( p.lU| ) and ( |111| ) if the functions hi{x'^, z) are known as solutions 
of the equations ( |109| ). 



An elongated rotation ellipsoid hypersurface is given by the formula |2T 



+ -u^ z^ 
— H 



P 



(112) 



where a > 1 and p is similar to the radial coordinate in the spherical symmetric 
case. 

The space 3D coordinate system is defined 

X = psinhMsinvcos(yO, y = psinhusinv sirup, 
z = p cosh M cos f, 

where a = coshw, (0<m<oo, < v < it, < f < 27r). The hypersurface metric 
is 



gvv = p 



(^sinh^ u + sin^ 



dipip = P^ sinh^ u sin^ v. 
Let us introduce a d-metric 



(113) 



6s = gi{u,v)du + dv + h^{u,v,ip) {6t) + h4{u,v,ip) {dipY 



(114) 



where St and Sip are N-elongated differentials. 

As a particular solution ( |108| ) for the h-metric, considering Tg = Tf, we choose 
the coefficient 

gi{u,v) = cos^ V. (115) 
The hs{u,v,ip) = h^{u,v,p{u,v,ip)) is considered as 



h3{u,v,p) 



1 



1-^ 



sinh u + sin v 



l + l£ 

4p. 



6 • 



(116) 



In order to define the coefficient solving the Einstein equations, for simplicity, 
with a diagonal energy-momentum d-tensor for vanishing pressure, we must solve 
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the equation ( |1U9D which transforms into a hnear equation if Ti = 0. In our case 
s {u, V, if) = (3~^ {u, V, if) , where (3 = {dh4^/d{p) /h4^, must be a solution of 



ds ^ dln^/lhsl 



dip dip 2 
After two integrations (see [|T^) the general solution for hi{u,v,p), is 

F{u, V, z) dz 



hiiu^ V, p) = 04 (m, v) exp 



;ii7) 



where 



F{u,v,z) = {^J\hs{u,v,z)\[sno) {u,v) + - / ^J\hs{u,v, z)\dz])~'^, ( 



:il8) 



zo{u,v) 



si(o) (m, v) and zq (m, v) are some functions of necessary smooth class. We note that 
if we put /i4 = a4{u, v) the equations (|109|) are satisfied for every = h^i^u, v, ip). 

Every d-metric (|114|) with coefficients of type (|115|) , (|116|) and ( |117| ) solves the 
Einstein equations (|108|) -( pTl| ) with the diagonal momentum d-tensor 

= diag [0, 0, —e = —mo, 0] , 

when Tg = 2Kmo; we set the light constant c = 1. If we choose 



04 (m, v) 



sinh^ u sin^ v 
sinh^ u + sin^ v 



our solution is conformally equivalent (if not considering the time-time component) 
to the hypersurface metric ( p.l3| ). The condition of vanishing of the coefficient ( |116| ) 
parametrizes the rotation ellipsoid for the horizon 



— H 

- 1 ^2 



(119) 



where the radial coordinate is redefined via relation r = p (l + ^) . After multipli- 
cation on the conformal factor 



^sinh^ u + sin^ 



1 + ^ 
4p 



approximating gi{u,v) = sin^ v ~ 0, in the limit of locally isotropic spherical sym- 
metry, 



the d-metric ( |114| ) reduces to 



ds' 



1 + ^ 
4p 



X +y + z = r 



(dx^ + dy"^ + dz' 



4p 



1 + ^ 
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which is just the Schwazschild solution with the redefined radial coordinate when 
the space component becomes conformally Euclidean. 

So, the d-metric ( |114| ), the coefficients of N-connection being solutions of (|110|) 
and ( |111| ), describe a static 4D solution of the Einstein equations when instead of a 
spherical symmetric horizon one considers a locally anisotropic deformation to the 
hypersurface of rotation elongated ellipsoid. 
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